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Abstract 

The femtosecond response of NO-doped rare gas matrices is studied within a stochastic Langevin 
theoretical framework. As is shown, a simple damped harmonic oscillator model can describe prop- 
erly the absorption and emission line shapes associated with the NO (A^T,^ < — > X^II) electronic 
transitions inside these media as well as the matrix first-solvation shell response in a process with 
two timescales, finding a fairly good agreement with available experimental data. This approach 
thus constitutes an alternative and complementary way to analyze the structural relaxation dy- 
namics of systems in liquids and solids, leading to a better understanding of the underlying physics. 
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Many molecular systems undergo local structural changes when they are photoexcited 
to electronic states with strong electron-phonon couplings. These structural changes may 
lead to an important reorganization of the surrounding medium [l, 2|- In this regard, due 
to their simple structural and well-known physical and thermodynamical properties, pure 
and doped rare gas matrices constitute an ideal benchmark to analyze and understand 
the fundamental mechanisms of configurational rearrangements in more complex molecular 
systems 

mm 

. Among them, NO-doped raregas matrices result of particular interest 
due to the many conventional (steady-state) 

mm 

, [lO|, lUl and femtosecond pump-probe 
spectroscopy experimental studies available for these systems, which 
reveal the medium dynamical response when the NO molecule is photoexcited to one of its 
lowest Rydberg states. 

From previous studies in the literature , the dynamics associated with the 

NO electronic transition between the ground state (X) and one of its lowest Rydberg excited 
states (e.g.. A) can be described as a four-step process , [lo| within the configuration 
coordinate model in the harmonic approximation (CCM-HA), which allows us to relate exper- 
imental spectroscopic information to the potential energy surfaces of the different electronic 
states of the impurity. The four-step process is essentially as follows. In a first step, the 
absorption of a photon excites the only unpaired electron of the NO molecule to a Rydberg 
state. This provokes an overlapping between the NO Rydberg orbital and the matrix shells 
constituted by the neighboring rare gas (Rg) atoms due to the larger effective size of the 
new NO excited state with respect to its ground state. As a consequence, the atoms of the 
first solvation shell undergo a strong repulsion, which displaces them from their equilibrium 
positions and, therefore, leads to an increase of the NO-matrix bond lengths. Because the 
process is energetically unfavorable, in a second step, the lattice relaxes around the excited 
impurity, with the Rg-atoms reaching a new equilibrium position. During this process, most 
of the energy in excess is dissipated through the matrix, but a small fraction is invested into 
an expanssion of the surroundings beyond the first solvation shell. This latter amount of 
energy is going to contribute substantially to the restoring force in the next step. This ba- 
sic mechanism consisting of an oscillating radial expansion/compression of the cage around 
the NO molecule, from one equilibrium cage radius to another one, is known as the elec- 
tronic bubble formation 

mm 

, and is also operative in rare gas liquids and clusters 
For instance, molecular dynamics simulations of this bubble formation 
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have been carried out 



21 



22 



23| to interpret the absorption and fluorescence spectra of 



Xe-dopped Ar^ , clusters upon photoexcitation of the Xe dopant to low-n Rydberg states 



171, m \m 

|20|. In the case of photoexcited impurities trapped in rare gases and van der 
Waals sohds, an extensive hterature can also be found lia, ll3|, l2J,l25|. In 

a third step, Rydberg fluorescence with the emission of a photon may occur, which (fourth 
step) brings the system back to its initial ground state conflguration by relaxation, leading 
to the equilibrium cage radius associated with the NO ground state. 

As infers from the above described absorption/emission process, the flrst and third steps 
are related to energetic aspects of the transition (steady-state spectroscopy), while the sec- 
ond and fourth steps are associated with the solvent relaxation dynamics (femtosecond 
pump-probe spectroscopy). In general, all rare gas matrices respond to photoexcitation 
by undergoing an ultrafast expansion of the flrst solvation shell during the flrst 100-200 fs 
(something similar also happens in the case of the emission process), where almost all the 
solute-solvent interaction energy becomes kinetic. During this expansion, the cage radius 
reaches a maximum above ri'^^ , the equilibrium cage radius for the NO molecule in its Ryd- 
berg excited state, and then gets back below it (but always above rlq , the equilibrium cage 
radius for the NO molecule in its ground state). The time at which the cage radius reaches 
this flrst minimum is the so-called first recurrence time t [ll, Q, Q, 1^. Due to the energy 
relaxation process, for any time t > t, the cage radius fluctuates with decreasing amplitude 
around ri'^\ which is flnally reached asymptotically (like a damped motion). 

The ultrafast lattice structural dynamics mediated by NO excitation has been studied in 



the last years by means of classical molecular dynamics (CMD) 26 



28 



29|, obtaining 



a reasonable agreement with the experimental data. Apart from numerical issues (e.g., 
introduction of periodic boundary conditions, cutoff radii or smoothing functions to avoid 
'border effects'), an important source for deviations between CMD-based simulations and 
experiment j^, [lo| is the unavailability of accurate interatomic potential surfaces describing 
the NO(74)-Rg interaction, where Rg = Ne, Ar, Kr or Xe. Therefore, the standard procedure 
to obtain the potential surface is by considering pairwise potential models and then carrying 
out simulations where the parameters deflning such models are varied until an optimal flt 
to the experimental data is found. Non-equilibrium and dissipative properties involved in 
this kind of processes (e.g., dissipation, relaxation or friction) can also be studied by means 
of CMD. Nevertheless, besides CMD simulations, these properties can also be extracted 
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alternatively by means of a stochastic approach, which renders a complementary view of the 
dynamics underlying the corresponding experimental findings. 

Here, we propose a fully stochastic classical Langevin model (CLM) to tackle the prob- 
lem of the analysis and interpretation of the lattice relaxation dynamics of NO-doped rare 



L ayn 



n 



gas matrices. More specifically, making use of the CCM-HA UMM, 

where the pairwise 

impurity-lattice and lattice-lattice interactions are substituted by a global cage harmonic 
potential, the structural relaxation of the first solvation shell of the solid matrix is described 
in terms of a simple damped harmonic oscillator (DHO) model, where the energy transfer 
to/from the remaining shells of the rare gas matrix (heat bath) is accounted for by a damping 



rate 



30 



3l| . In this way, the DHO-CLM used here avoids the problem of considering ex- 



plicitly particular many-body potential energy surfaces, this resulting in a gain of analytical 
handiness and, therefore, of physical insight. Moreover, it also entails a considerable save of 
computational time, something very important when one is interested in the global under- 
standing and fitting procedure of the experimental data rather than in a detailed dynamics 



accounting for them. In this regard, C 



study surface diffusion [32 
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35, 



iM-based approaches have been successfully used to 



36(1 . It is worth stressing that, though the Langevin 

'ormally derived from 



equation employed in CLMs is seemingly phenoraen olog ical, it can be 

39, 3: 



the so-called Caldeira-Leggett Hami 
the theory of open quantum systems 



toman 
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43l |. very well known in 



3l| . The Caldeira-Leggett Hamiltonian, introduced 



in the dynamics of open classical/quantum systems in 1983, shows that the effects of the 
detailed dynamics of a given bath or thermal reservoir (solid, liquid, etc.) acting on a system 
can be replaced by an infinite sum of harmonic oscillators. The damping acting over the sys- 
tem can be expressed, in general, as a time-dependent function involving the corresponding 
harmonic frequencies and the system-bath coupling constants. When the damping function 
becomes Ohmic (i.e., it is described by a constant spectral function of the frequency), the 
generalized Langevin equation (arising from the Caldeira-Leggett Hamiltonian) describing 
the system evolution reduces to a standard Langevin equation. 

In the problem of NO-doped rare gas matrices, the electron-phonon coupling during the 
electron excitation has been estimated from the experiment [7] and is directly related to 
the displacements of the equilibrium configuration coordinate in the ground and excited 
states. On the other hand, the lattice relaxation (first solvation shell) can be described by 
a new damping rate which gives us an idea of the energy loss to the Rg-atom matrix acting 
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as a heat bath via a sort of effective phonon-phonon couphng after excitation within the 
electronic bubble model. The dynamics of the cage radius involving the first solvation shell 
is coupled to the remaining solvation shells displaying a damped harmonic oscillation. As 
will be shown below, this damping rate introduces a new timescale. 

Thus, the NO absorption/emission spectrum is related to the response of the matrix first 
solvation shell (i.e., the shell of first nearest neighbors around the impurity). In order to 



describe the dynamical 



evolution of the first solvation shell, the CCM-HA is commonly used 



in the literature [3, S 1^, with the configuration coordinate being the matrix cage radius^ 
r. Within this model, r represents the mean distance between the NO impurity and each 
one of the atoms of the first solvation shell, and an effective harmonic potential (per mass 
unit), 

U{r) = \ujl{r-r,,)\ (1) 
stands for the impurity- matrix interaction. In ([1]), ujq is the effective harmonic frequency 



and Teq is the equilibrium cage radius of the first solvation shell 1271 . bsl . |29 |: 

fx) 

and Teq will denote the equilibrium cage radii after absorption and emission, respectively. 
Estimates of the harmonic frequency and the cage radius increment, A = ri^'' — have 
been obtained in the literature for different matrices from a direct analysis of experimental 
spectroscopic data by using either the method of moments f^j or the semiclassical projection 
method [o], Q. The so-called Huang- Rhys factors [7] suggest a decreasing electron-phonon 
coupling with the rare gas mass forming the solid matrix. This can be understood by taking 
into account the fact that, as the Rg atom mass increases, the characteristic NO-Rg bond 
lengths, and also the Rg-Rg ones, become larger and, therefore, the NO Rydberg orbital can 
be better accommodated within the space inside the first solvation shell. The values of the 
effective harmonic frequency in the excited state, as well as the cage radius increment 
found for NO-Rg matrices, with Rg = Ne, Ar, Kr and Xe, are given in Table [H As can be 
noticed, while A decreases monotonically as the NO-Rg and Rg-Rg bond length increase for 
ujo, the variation with the mass follows the opposite trend. Notice that an increase of ujq with 
the mass results counterintuitive, for one would be tempted to think that heavier particles 
will oscillate slower than lighter ones. As will be shown, in our model the ujq frequencies, 
as well as the corresponding recurrence times, will be modified due to the presence of the 
damping rate. 

After assuming the CCM-HA, the effects of the energy transfer process between the first 



(X) 

TABLE I: Effective liarmonic frequencies associated witli the NO in the ground {lOq ) and the 
Rydberg excited {^^q^^) states, and increment of the NO-Rg cage radius (A) for Rg = Ne, Ar, Kr 
and Xe, within the configuration coordinate model in the harmonic approximation. For Ne, Ar 
and Xe, the value of these magnitudes have been taken from Ref. |lo|, while for Kr they are from 
Ref. [7|. 

Ne Ar Kr Xe 





(A) 



(ps-1) 1.26 1.23 1.35 1.53 

(ps-^) 1.74 1.41 1.45 1.59 

A (A) 0.59 0.41 0.18 0.14 



solvation shell and the remaining shells from the solid matrix can be described in terms of 
the cage radius by means of the Langevin equation associated with a DHO, as 



f(t) = -7r(t)+F[r(t)]+eW- 



(2) 



In this equation, F = —VU is the deterministic force (per mass unit) acting over the 
first solvation shell. The fluctuating force, ^(t), is described by a Gaussian white noise, 
thus satisfying {^(t)) = and {C{t)C{t')) = 2fi'ykBT6(t — t'), where /i is the effective mass 
associated with the first solvation shell, i.e., its value is 12 times the Rg-atom mass 
lOl], T is the temperature and ks the Boltzmann constant. Furthermore, 7 is a constant 
damping rate related to the fluctuating force ^(t) through the fluctuation- dissipation theorem 



30, l3l|- As can be noticed, Eq. ([2]) is a standard Langevin equation, which is used instead 



of its generalized version because the relaxation time at which the system (the cage radius) 
approaches equilibrium with the heat bath is much longer than the correlation time of the 
bath fluctuations {Markovian approximation) . Phenomenological values of 7 have been used 
in the literature in stochastic Langevin theoretical treatments to describe solute-solvent 
energy exchange in liquids and solids 4J, |45|, |46|, |47| , though it is more commonly related 



to liquids and gases. 

In our simulations, 10,000 stochastic trajectories have been run for each electronic state up 



to 5 ps, considering time steps of 



CMD simulations 



26 



271, 



28 



fs. Initial conditions are assigned like in the corresponding 
29[. The initial positions are randomly generated according 



to a uniform distribution centered around rif^ or Ve^' , depending on whether emission or 



TABLE II: Optimal values of 7 obtained from the DHO-CLM. This damping rate governs the 
relaxation dynamics in the ground state after emission (7x) or the excited state after absorption 
(7a). 





Ne 


Ar 


Kr 


Xe 


ix (ps ^) 


0.51 


0.57 


0.67 


1.15 


lA (PS~^) 


0.51 


0.61 


1.05 


1.55 



absorption are being studied, respectively. Regarding the choice of the initial velocities, 
they follow a standard Boltzman distribution at a temperature T. The experiments are 
carried out at T = 4 K, where the nuclear dynamics of the matrix atoms is dominated by 
the zero-point motion. Hence, in order to account for this quantum effect properly in CMD 
simulations, a higher effective temperature depending on the Debye frequency of the solid 



has to be considered 



48l | (e.g., for a solid Ar matrix, the effective temperature corresponding 



to T = 4 K is about 49 K). Thus, the initial velocities have to be re-scaled accordingly (i.e., 
taking into account this effective temperature) before starting the simulation, which we have 



carried out by using the velo city- Ver let algorithm |49| to propagate in time the stochastic 
trajectories. Then, the absorption and emission bands are calculated by constructing his- 
tograms with the energy differences between the initial (ground state for absorption, excited 
state for emission) and the final state (excited state for absorption, ground state for emis- 
sion). In particular, we have considered the electronic transition between the NO ground 
state and its first lowest Rydberg state, i.e., A^S^ < — > X^II. The medium response to 
the electronic transition is induced by an instantaneous switching from the ground state 
potential to the excited one (making use of the Franck-Condon principle). Now, in order 
to reproduce the experimental bandwidths, 7 is chosen as the only fitting parameter (re- 
member that in the homologous CMD simulations, the fitting parameters usually appear in 
the potential energy surface models considered), which is varied until the band widths are 
fairly reproduced. The optimal 7 values obtained from the DHO-CLM are given in Table HTl 
These values deviate less than 10% from the rates obtained through CMD after computing 
the approximated time exponential decay of the total energy autocorrelation function for 
the impurity, 

CE{t) = m)m-t^))-{E\t)). (3) 
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From this expression one can extract the rate of energy transfer between the first solvation 
shell and the surrounding bath during the corresponding relaxation process and, therefore, 
of the cage radius damping towards its equilibrium value. Fortunately, the band widths 
are not very sensitive to the 7 rates allowing us a certain degree of flexibility in the fitting 
procedure. As can be seen, 7 increases from Ne to Xe matrices and is always larger in 
the excited than in the ground state. Physically, this indicates an important energy loss 
between the first solvation shell and the medium or heat bath via an effective phonon- 
phonon coupling, which increases with the mass of the species constituting the solid matrix 
and, therefore, its rigidity. Moreover, this energy loss is also larger when the impurity is 
photoexcited. In the experiment and also in CMD simulations, a rapid relaxation of the 



26 



28 



29[ for Ar, Kr and Xe 



photoexcited NO molecules to the medium is observed 
matrices, while for Ne matrices the effect is slower and "collective" (see below). Finally, as 
Eq. ([2]) describes the dynamics of the first solvation shell, 7"^ will also provide a timescale 
for the solvent relaxation dynamics, which ranges from ~650 fs for NO (A) in Xe to ~2,000 fs 
in Ne. 

In Fig. [U absorption (right/blue) and emission (left/red) bands obtained from the DHO- 
CLM simulations are plotted for matrices constituted by different Rg-atoms and compared 
with the corresponding experimental data (dots). These steady-state absorption and emis- 
sion bands represent the initial and final stages of the photoinduced structural rearrangement 
process around the excited NO molecules. To compare with the effects produced by the pres- 
ence of the matrix, the green horizontal dashed line around 5.5 eV indicates the position of 
the NO absorption/emission electronic transition in the gas phase. As can be noticed, the 
line shapes of the bands and their Stokes shifts are in good agreement with the experimental 



data, the main differences being lower than 2% in al 
can be found, however, in the CMD simulations 26, 



cases. Larger numerical discrepancies 



27l |. For example, for Ar matrices, the 



DHO-CLM reproduces quite accurately the experimental value of the Stokes shift, 0.570 eV 







while the CMD calculations provide a shift of 0.700 eV 26|, l27|]. In the case of Ne, the 
shift obtained from CMD is not so good in comparison with the experiment and depends 
strongly on the parameters of the pairwise potentials considered. On the other hand, al- 
though for Kr and Xe matrices both the DHO-CLM and CMD [29] simulations reproduce 
fairly well the experimental value of the Stokes shift [3], the corresponding line shapes have 
not been calculated from the latter. 
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FIG. 1: (Color online.) Absorption (right/blue) and emission (left/red) line shapes associated 
with the NO <-> X^H) electronic transition in rare gas solids. In all cases, experimental 

data are represented by dots, while the line shapes obtained with the DHO-CLM are denoted 
by solid lines. The green horizontal dashed line around 5.5 eV denotes the position of the same 
absorption/emission process, but for the NO molecule in gas phase. 

The time evolution of the cage radius after the absorption is plotted in Fig.[2]^a), where the 
cage radius is expressed in terms of its effective increase, r' = r — req . In this way, at t = 0, 
r' = means that the cage radius starts at r = Veq , on the contrary, asymptotically, in the 
limit t oo, r' will approach the cage radius increment A, since the cage radius approaches 
the equilibrium value, rif\ Physically, this curve describes a process with two different 
timescales. The first part of the process consists of the initial cage radius expansion (from 
t = to the first maximum) and indicates the matrix response after a photon is absorbed by 
NO and its unpaired electron is promoted to the Rydberg state, thus increasing the effective 
size of the impurity (and so the size of the bubble). As seen in Fig. EJ^a), this expansion 
(solid line) is relatively fast and mimicking the CMD results (blue dotted line) for about the 
first 150 fs in all cases. In particular, the maximum expansion velocity (slope) reached by 
the first solvation shell is 471 m/s (139 fs) for Ne-matrices, 308 m/s (171 fs) for Ar-matrices, 
167 m/s (163 fs) for Kr-matrices, and 114 m/s (141 fs) for Xe-matrices (between brackets, 
the time at which this velocity is reached), which are similar to those obtained from CMD 
[l^ . However, because the electron-phonon coupling decreases with the Rg-atom mass, as 
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derived from the method of moments in Ref. [3], not only a smaller increase of the cage 
radius from Ne to Xe matrices is noticeable, but also in the initial slope of its evolution, as 
seen in Fig. I21^b). The second part of the process describes the relaxation of the cage radius 
to its equilibrium value after the initial boosting led by the (ultrafast) expansion of the NO 
orbital. This is the stage that follows the first maximum in the curves displayed in Fig. [2](a). 
As can be noticed, the damping of the oscillations around the corresponding equilibrium cage 
radius increases from Ne to Xe matrices, this being indicated by an also increasing value 
of 7 (which we have obtained from the fitting to the experimental spectroscopic measures, 
as indicated above). According to the DHO model, the corresponding oscillations displayed 
in Fig. [2] and their envelopes follow an expression of the type e'^^l"^ cos(a;it + 5), where the 



renormalized harmonic frequencies are given 



30| by 



and the phase shift by 5 = tan~^(7/2c(}). From Ne to Xe matrices, these renormalized 
frequencies are 1.72, 1.38, 1.35, and 1.39 ps~^, respectively. These values follow the same 
tendency observed in the experiment 7|] by means of an analysis carried out by the method 
of moments (see Table [T]). In principle, a larger damping implies a lower oscillation, which 
stresses with the fact that the Rg-atom mass increases. In other words, intuitively, one 
should expect that the heavier the Rg-atom mass, the slower the oscillation and, therefore, 
the larger the recurrence times (which are discussed below). In order to understand this 
counterintuitive behavior, remember that bond lengths increase with the Rg-atom mass 
and, therefore, one could expect a matrix response more typical of rigid solids for larger Rg 
atoms, and more "floppy", i.e., involving more neighboring shells, for smaller Rg atoms. If 
so, as happens with sound waves, the transmission of energy through the material is more 
effective (faster) for solids than for liquid which, within the context of our model, means 
greater energy losses to the first-neighbors shells as the Rg-atom mass increases. 

As shown in Fig. [2]^a), for Ne matrices the DHO-CLM renders a value for the first 



recurrence time o: 
Texp ^ 1.1-1.4 ps 



7pgn_cLM ~ 564 fs, which does not agree with the experimental value 
16j. In the case of CMD, in order to obtain a good estimate of this 



15|, 



time, an artificially shortened Born-Mayer potential model was considered, which provides 



TcMD ~ 1.3 ps [l5|, l28|. To understand this discrepancy, as mentioned above, note that 



the NO(X)-Ne and Ne-Ne bond lengths are smaller than the NO(74)-Ne one, which covers 
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FIG. 2: (Color online.) Cage radius dynamics, a) Blue dotted lines have been obtained from 



CMD simulations 



271, 



28 



2^, while black solid lines represent the results from the DHO-CLM. b) 
Comparison of the four cases considered in part (a): Ne-matrix (black solid line), Ar- matrix (blue 
solid line with circles), Kr-matrix (green dashed line) and Xe-matrix (red dotted line). 



about four solvation shells 



15 



28| . Thus, although the initial expansion of the cage radius 
is well described by the DHO-CLM because it is only affected by the first solvation shell, 
immediately after such an expansion the size of the excited NO wave function overlaps the 
matrix atoms of upper shells, which are not considered in the DHO-CLM. The contribution 
of these upper shells would lead to observe the first recurrence at longer times due to a slower 
relaxation (i.e., a "floppy" -like matrix behavior under a local structural perturbation). The 
recurrence time, therefore, would not be ruled by one shell but by four shells at a time, as 



shown through CMD [15|, |28|]. A simple way to deal with this effect consists of including 
an anharmonic contribution in the effective potential acting over the solvation shell in the 
involved excited state. This procedure improves the value of the recurrence time to be 
tdho-clm ~ 1-09 ps and leaves unchanged the corresponding line shapes. Obviously, this is 
the simplest and fastest way of considering the inclusion of more external shells. Other more 
elaborated schemes can be envisaged such as, for example, using a four-variable model with 
four coupled Langevin equations, each one describing the evolution of the radius associated 
with each shell. In the case of Ar matrices, we find a recurrence time greater than in the 
corresponding CMD calculation [see Fig. [2](a)], tcmd ~ 500 fs [l^ versusTpno-CLM ~ 711 fs, 
in better agreement with the experimental value Texp ~ 600-800 fs Q, [iGj. Contrary to the 
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case of Ne, here the N0(y4)-Ar bond length is such that the expansion of the NO wave 
function up to t ~ r is going to affect essentially the first solvation shell and, therefore, the 
dynamics will be well reproduced. Note that, as /i increases (from Ne to Xe), the amplitude 
of the expansion curve of the cage radius decreases as well as the expansion velocity of the 
solvent shell. However, the first recurrence time is not an increasing monotomic function of 
the mass, as observes in Fig. [2](b). With respect to Kr matrices, we have found an acceptable 
similarity in the response between DHO-CLM and CMD [note in Fig. [2](a) that tdho-clm ~ 
tcmd]- For Xe matrices, both simulations provide a similar behavior pattern qualitatively 
(regarding the amplitude of the oscillations) though some quantitative discrepancies can be 
found. To our knowledge, there are no reported femtosecond experiments for the bubble 
dynamics in Kr and Xe matrices and, therefore, the corresponding curves (as also happens 
with those obtained from CMD) have to be considered as predictive. 

Summarizing, we have shown that the experimental absorption/emission bands associated 
with the A^Ti^ ^ X^U electronic transition of NO in rare gas matrices can be fairly well 
reproduced by means of a simple damped harmonic oscillator model, which describes the 
dynamics of the first solvation shell (i.e., the response of the matrix to the NO excitation) in a 
process with two different timescales. The first part of the process (fast expansion of the first 
solvation shell) is directly related to the electron-phonon coupling which has measured in the 
experiment. In the second part of the process, the introduction of a new damping rate via an 
effective phonon-phonon coupling after excitation in order to describe the lattice relaxation 
allows us to obtain very reliable band structures by only using data from spectroscopic 
measurements and a single fitting parameter (7). As shown, the expansion dynamics of 
the cage radius is well described for about the first 200 fs as well as the first recurrence 
time. Our DHO-CLM should improve with the mass of the matrix atoms. Note that the 
calculation of r from the model presented here marks the time at which a well-defined first 
recurrence takes place. Nonetheless, the fact that only the first solvation shell is involved 
in the dynamics described restricts the applicability of the model when 'collective' shell 
behaviors are present, as in the case of Ne matrices. In such cases, it has been shown that the 
addition of anharmonic contributions to the effective harmonic potential used, which would 
account for the eventual distortions undergone by the first solvation shell, render better 
recurrence times in comparison with the experimental ones. For Ar matrices, tdho-clm 
agrees fairly well with the experimental value without any need to consider anharmonicities. 
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The oscillations of the cage radius as well as the recurrence times are easily explained 
within the DHO model through a renormalization of the harmonic frequencies due to the 
presence of the damping rate. This information could be therefore used to design new 
femtosecond pump-probe experiments. The fact that all of these features arc reasonably 
well reproduced by simple, stochastic Langevin simulations opens up a new theoretical way 
to explore, describe and understand the physics underlying rearrangement processes in solids, 
stressing the insightful value of the model. 
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